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O Abstract 

gn! We assume that the Pauli exclusion principle is violated for neutrinos, and thus, 

(— I ' neutrinos obey at least partly the Bose-Einstein statistics. The parameter sin 2 x 

is introduced that characterizes the bosonic (symmetric) fraction of the neutrino 
, wave function. Consequences of the violation of the exclusion principle for the two- 

neutrino double beta decays (2z^/3/3-decays) are considered. This violation strongly 
changes the rates of the decays and modifies the energy and angular distributions of 
. the emitted electrons. Pure bosonic neutrinos are excluded by the present data. In 

the case of partly bosonic (or mixed-statistics) neutrinos the analysis of the existing 
data allows to put the conservative upper bound sin 2 \ < 0-6- The sensitivity of 
future measurements of the 2i//3/3-decay to sin 2 \ is evaluated. 

o \ 

! 1 Introduction 

X' 

Does neutrino respect the exclusion principle of it's inventor? In this paper we assume 
that Pauli exclusion principle is violated for neutrinos and therefore neutrinos obey (at 
least partly) the Bose-Einstein statistics. Possible violation of the exclusion principle was 
discussed in a series of papers p] though no satisfactory and consistent mechanism of the 
violation has been proposed so far. The assumption of violation of the Pauli exclusion 
principle leads to a number of fundamental problems which include loss of a positive defi- 
niteness of energy, violation of the CPT invariance, and possibly, of the Lorentz invariance 
as well as of the unitarity of S-matrix. (For a critical review see ref. [2].) Experimental 
searches of the effects of the Pauli principle violation for electrons [3] and nucleons [I] have 
given negative results, leading to extremely strong bounds on the magnitude of violation. 



^1" 
^1" 



1 



It may happen however that due to unique properties of neutrinos (neutrality, small- 
ness of mass associated to some high mass scales), a violation of the Pauli principle in the 
neutrino sector is much stronger than in other particle sectors. Therefore one may expect 
that effects of its violation can be first seen in neutrino physics. 

A possibility of the Bose statistics for neutrinos has been first considered in ref. [5] 
where its effects on the Big Bang Nucleosynthesis (BBN) have been studied. According 
to [5] the change of neutrino statistics from pure fermionic to pure bosonic diminishes the 
primordial 4 He abundance by ~ 4%. 

The idea of bosonic neutrinos has been proposed independently in ref. [6], where 
cosmological and astrophysical consequences of this hypothesis have been studied. Bosonic 
neutrinos might form a cosmological Bose condensate which could account for all (or a part 
of) the dark matter in the universe. "Wrong" statistics of neutrinos modifies the BBN, 
leading to the effective number of neutrino species smaller than three. The conclusion 
in [6] agrees qualitatively with results of [5] though quantitatively a smaller decrease of 
N v is found [7]. 

As far as the astrophysical consequences are concerned, dynamics of the supernova 
collapse would be influenced and spectra of the supernova neutrinos may change [6j [8] . 
The presence of neutrino condensate would enhance contributions of the Z-bursts to the 
flux of the UHE cosmic rays and lead to substantial refraction effects for neutrinos from 
remote sources [6]. 

We assume that the Pauli principle is violated substantially for neutrinos, while the 
violation is negligible for other particles. In particular, for electrons we will assume the 
usual Fermi-Dirac (FD) statistics. How to reconcile this pattern of the violation with the 
fact that in the standard model the left-handed neutrino and electron belong to the same 
doublet? The answer may be connected to the fact that neutrinos are the only known 
neutral leptons and thus they can have substantially different properties from those of the 
charged leptons. In particular, neutrinos can be the Majorana particles and violate lepton 
number conservation. The difference between charged leptons and neutrinos should be 
related to breaking of the electro-weak (EW) symmetry, and it can originate from some 
high mass scale of nature. One may consider scenario where violation of the Pauli prin- 
ciple occurs in a hidden sector of theory related to the Planck scale physics, or strings 
physics. It could be mediated by some singlets of the Standard model - (heavy) neutral 
fermions which mix with neutrinos when the EW symmetry is broken. Since only neutri- 
nos can mix with the singlets, effects of the Pauli principle violation would show up first 
in the neutrino sector and then communicate to other particles. In this way a small or 
partial violation of the relation between spin and statistics might occur. A violation of 
the spin-statistics theorem for other particles can be suppressed by an additional power 
of a small parameter relevant for the violation in the neutrino sector and due to weak 
coupling of neutrino to other particle sector. 

A violation of the Pauli principle for neutrinos should show up in the elementary 
processes where identical neutrinos are involved. A realistic process for this test is the 
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two-neutrino double beta decay (2z//?/3-decay), 

A -> A' + 1v + 2e~ (1) 

(or similar with neutrinos and positrons). It was shown in [6] that the probability of the 
decay as well as the energy spectrum and angular distribution of electrons should be af- 
fected. Qualitative conclusions were that the pure bosonic neutrino is excluded, whereas 
large fraction of the bosonic component in a neutrino state is still allowed by the present 
data. In this connection, a possibility of partly bosonic (mixed-statistics) neutrinos can 
be considered. 

In this paper we perform a detailed study of the effects of bosonic neutrinos on the 
double beta decay. In sect. 2 we consider the general case of partly bosonic neutrinos. 
We introduce a phenomenological parameter sin 2 \ which describes the fraction of bosonic 
neutrinos in such a way that a smooth change of sin 2 \ from to 1 transforms fermionic 
neutrinos into bosonic ones. So, in general, neutrinos may possess a kind of mixed or more 
general statistics than Bose or Fermi ones [HI EI]. In sect. 3 we present an analytic study 
of the double beta decay probabilities. The exact expressions for the 2z//3/3-decay rates to 
ground and excited + and 2 + states with corresponding nuclear matrix elements (NME's) 
are given in sect. 4. The results of numerical calculations of the total rates and various 
distributions for the 2i//?/9-decays of 76 Ge and 100 Mo are presented in sect 5. In sect. 6. 
we obtain the bounds on sin 2 \ from the existing data and evaluate the sensitivities of 
future double beta decay experiments. Discussion and conclusions are given in sect. 6. 

2 The 2z//3/3-decay for bosonic and partly bosonic neu- 
trinos 

In the case of mixed statistics the operator of neutrino state can be written as 

\v) = a + \0) = c s f + \0) + sib + \0) = c s \f) + s s \b) (2) 

where |/) and |6) are respectively one particle fermionic and bosonic states. The normal- 
ization of \v) implies c 2 + s 2 = 1 (c§ = cos 5 and sg = sin 5). f (f + ) and b (b + ) denote 
fermionic, and bosonic annihilation (creation) operators. 

To develop a formalism for description of identical neutrinos one needs to specify 
commutation/anti-commutation relations. We assume that they have the following form: 

fb = e%f 1 f + b + = e% + f + , fb + = e-% + f, f + b = e~%f + , (3) 

where <fi is an arbitrary phase. Then the two-neutrino state can be defined as 

Ifci.Afc) = a+o+|0). (4) 

For the pure bosonic neutrino one cannot introduce the Majorana mass term. So, the 
neutrinoless double beta decay should be absent. In the case of partly bosonic neutrino, 
the neutrino mass would appear due to its fermionic component. This means that the 
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kinematical mass measured, e.g. in the tritium beta decay, would not be the same as the 
mass found from the neutrinoless beta decay. Such a situation, however, can be realized 
in the case of the usual fermionic neutrinos too. 

The amplitude of the decay of nucleus A — > 2v + 2e + A' can be written as 

A 2 p = (e(p e i),e(p e2 ), F(p„i), V(p v2 ), A'\ J d 4 x 1 d 4 x 2 '4> l/ (x 1 )?p l ,(x2)M( y x 1 , x 2 )\A). (5) 

After making the necessary commutation, according to eq. ([H]), we obtain 

A w = A f [4 + c 2 s 2 (1 - cos 0)] + A b [sj + c 2 s 2 (1 + cos 0)] , (6) 

where Af and Af, are respectively fermionic (antisymmetric) and bosonic (symmetric) 
parts of two antineutrino emission. The amplitude can be parametrized as 

A 2f3 = cos 2 xA f + sin 2 x A b , (7) 

where cos 2 x — c s + c 2 s 2 (1 — cos 4>) and sin 2 x = sj + c 2 s 2 (1 + cos (j>). 

After integration over the neutrino phase space an interference between fermionic 
Af and bosonic A b parts of the amplitude A 2 p vanishes because the fermionic part is 
antisymmetric with respect to neutrino interchange, while bosonic is symmetric. The 
probability of the 2z/ / 9 / 9-decay is equal to: 

W tot = C os 4 X Wf + sm A xW bl (8) 

where Wf tb are proportional to |y4y fe| 2 . The expressions for Wf jb will be given in the next 
section. 

Qualitative features of the (3(3— decay in the presence of the bosonic or partly bosonic 
neutrinos can be understood using the following consideration. Essentially, the effect of 
neutrino "bosonization" is that two contributions to the amplitude of the decay from 
diagrams with permuted neutrino momenta p v \ <-> p u2 should have relative plus sign 
instead of minus in the FD-case. 

The decay probability, W b , is proportional to the bilinear combinations of the type 
K h m K h n) K h m L h n , L b m L b n (see the next section), where 

Km = [Em — Ei + E e \ + Evi] 1 — [E m — Ei + E e2 + E v2 \~ x , 
L m = [E m ~ Ei + E e2 + E^Y 1 — [E m — Ei + E el + E v2 \~ x . (9) 

Here Ei is the energy of the initial nuclei, E m is the energy of the intermediate nuclei, 
E e j, and E U j are the energies of electrons and neutrinos respectively. The factors @ 
correspond to the propagators of the intermediate nucleus. The key difference between 
the bosonic and fermionic cases is the opposite signs of the two terms in the expressions 
Q). In the case of fermionic neutrinos they enter with the same signs (see, e.g. [1 lj ) : 

K f m = [E m - E t + E el + E^]- 1 + [E m -Ei + E e2 + E u2 ]~\ 
L f m = [E m -Ei + E e2 + E^}- 1 + [E m -Et + E el + E u2 ] " x . (10) 

(Remember that for electrons we assume the normal Fermi statistics.) The terms in (|9]) 
correspond to the amplitudes with permuted momenta of both neutrinos and electrons. 
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In the case of fermionic neutrinos such an interchange flips the sign twice (due to neutri- 
nos and electrons), so that the overall sigh turns out to be plus. In the case of bosonic 
neutrinos the permutation of electrons only changes the sign, and the overall sign is minus. 

Experimentally interesting are the 2z//3/3-decays to the ground states 0+ s and to excited 
states 0^ and 2f. The effect of bosonic neutrinos on the 2z//3/3-decay half-life is different 
for J 71 " = 2 + and J n = + . This can be understood qualitatively, approximating the 
combinations K h m and L b m for bosonic neutrinos by 



r b Ee2 — E e i + Ev2 ~ EjA. jb -Eel ~ E e2 + E v2 ~ E v \ 

^ (E m -E l + E /2y ' m " (E m -E l + E /2y ' 



and the corresponding combinations for the fermionic neutrinos by 

2 

Iff ~ tj ~ do) 

Here E /2 = (E e + E v ) is the average energy of the leptonic pair, E Q = Ei — Ef is the 
energy release in the decay, and Ef is the energy of the final nucleus. 

For the + — > + transitions an appearance of the differences of the electron and 
neutrino energies in the numerators of (TlTl) leads to substantial (1-3 orders of magnitude) 
suppression of the total probability. It also modifies the energy distributions of electrons. 

The effect of bosonic neutrinos on + — > 2 + transitions is opposite: The probabilities 
of transitions are proportional to the combinations [K h m — L^^K^ — L^), where 

( K b rb )r ^ 2{E e2 -E el ) 

{ m m) ~ (E m -E l + E /2r {l6) 

In the case of fermionic neutrinos the combination (K^ — L^J has an additional factor 
[E v2 — E v \)/(E m — Ei + Eq/2) and the suppression is stronger. Parametrically the prob- 
abilities of the + — > 2 + and + — > + transitions become of the same order for bosonic 
neutrinos. 

In the decay rates, the kinematical factors and L^ b are weighted with the corre- 
sponding nuclear matrix elements (NME's). Let us introduce the ratio 

ro{J ] = WfWr ( } 

of the decay probabilities to ground (J n = 0+ s ) and excited (J n = 0^, 2^) states in pure 
bosonic Wb(J 7T ) and pure fermionic cases Wf(J n ). In general, to find ro(J n ) one needs 
to calculate the NME for a given transition within an appropriate nuclear model. The 
situation is simplified for those nuclear systems, where the transition via solely the ground 
state of the intermediate nuclei m = 1 dominates [TJ1 EH E] . For those nuclei the single 
state dominance (SSD) approximation (hypothesis) can be used. In this case the NME's 
can be factored out in the rates and therefore cancel in the ratio ro(J n ). 

Let us consider the characteristics of the (3(3 decay to the ground and excited states J n 
in the mixed-statistic case of partly bosonic neutrinos. According to our considerations 



5 



the total decay probability and the normalized total differential rate can be written as 



where 



Wtat{r) = cos 4 x w f (.r) + sm 4 x w b (.r), (15) 

dW tot {r) = cos 4 xduJ f (Jn + sm 4 X ro(J 7T )duJ b (Jn 
[ ' W tot {J") cos 4 X + sin 4 xr (J-) ' 1 ' 

are the normalized distributions. Here dWf(J n ) and (iiy^J^) are the differential rates 
of the 2^/9/9-decay for the pure fermionic and bosonic neutrinos. In the case of single 
state dominance due to factorization, the normalized distributions do not depend on the 
uncertainties of the matrix elements [131 E]. I n general, the factorization does not occur 
and the uncertainties of nuclear matrix elements restrict substantially the sensitivity of 
the /3/3-decay to statistics of neutrinos. 



3 Rates and nuclear matrix elements 

For the cases of pure fermionic and bosonic neutrinos we outline the derivation of 2v/3f3- 
decay rates. The relevant nuclear matrix elements will be evaluated and discussed using 
the SSD and HSD (higher states dominance) hypothesis [131 Ej. 
The matrix element of the 2i//3/3-decay process takes the form 



i-i) 2 



<f\S®\i>= 

<e{p el ),e{p e2 ),V{p vl ),V{p v i),A!\T [H p {x 1 )H p (x 2 )] \A>dx 1 dx 2 , 



2 

where the weak /3-decay Hamiltonian is 

U^x) = ^ [e(x) 7 ^(l + 75>e(z)] -Ux) + h.c. (19) 

Here, J^{x) is the weak charged (nuclear) hadron current in the Heisenberg representation. 
The T-product of the two hadron currents can be written as 

T(H^x 1 )H^(x 2 )) = 
e(x 10 - x w )HP{x x )HP{x 2 ) + 6(x 20 - x 10 )H^x 2 )H^ Xl ). (20) 

In the derivation of the 2z//3/9-decay rate a number of conventional approximations have 
been used: i) Only the Si/ 2 wave states of the outgoing leptons are taken into account, ii) 
The contribution of the double Fermi matrix element to the decay rate is neglected as the 
initial and final nuclei belong to different isospin multiplets. iii) Only the leading order 
(1/rrip) Gamow- Teller operators in the non-relativistic reduction of the hadron current 
are retained. 
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For the differential 2i//3/3-decay rates to + ground state and 2 + excited state we obtain 

dW f , b (J + ) = a 2v F{Z h E el )F(Z f , E e2 ) M s fi d£l, (21) 

where a 2u = (C^A) 4 ^e/(64vr 7 ) and Gp = Gf cos 9 C (Gp is Fermi constant, 9 C is Cab- 
bibo angle). F(Zf,E e ) denotes the relativistic Coulomb factor and gA is the axial- vector 
coupling constant. The upper index / (b) stands for fermionic (bosonic) neutrinos. 
The phase space factor equals 



dQ = —E elPel E e2Pe2 E 2 X E 2 u2 5{E el + E e2 + E ul + E u2 + E f - x 
dE e i dE e2 dE v \ dE v2 dcos9. 



(22) 



Here, 9 is the angle between the outgoing electrons. (J n = + , 2 + ) consists of the 

products of nuclear matrix elements: 



with 



' 4 



i r /,6 , r f,b {2 _ 1 \rf> b _ rf> b \ 2 



Pel ■ Pe2 

E e \E e2 



'/,*>! 2 



1 1 Pel • Pe2 

3 E e iE e2 



v m j fc 

7=E < ^/iiE^-ii 1 - >< ^n£^A + > L - 



J+ 



(23) 



(24) 



Here, s = 1 for J = and s = 3 for J = 2. |0+ >, |0+ > (|2+ >) and |1+ > are, 
respectively, the states of the initial, final and intermediate nuclei with corresponding 
energies Ei, Ef and E m . The energy denominators and L(£ were introduced in Eqs. 
(ED and ([TOD. 



3.1 Higher states dominance 

The 2i//3/?-decay rates are usually evaluated in the approximation in which the sum of 
the two lepton energies in the denominator of the nuclear matrix element is replaced with 
their average value E /2 

E m -E t + E ej + E vk ~ E m — Ei + E /2 (25) 

(j,k = 1,2). The main purpose of this approximation is to factorize the lepton and 
nuclear parts in the calculation of the 2z//3/3-decay half-life. This approximation is justified 
if the transitions through the higher-lying states of the intermediate nucleus (at least few 
MeV above the ground state of (A,Z+1) nucleus) give the dominant contribution to the 
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2z/ / 3 / 3-decay amplitude. This assumption is called the higher states dominance (HSD) 
hypothesis. It is expected to be realized for A= 48, 76, 82, 130, 136 nuclear systems. 
Assuming the HSD hypothesis we obtain for fermionic neutrinos 



Ml = |A«(0+)f 



if I A/f( 3 ) (O+W 2 ~ Ee2) 2 (E u \ — E v2 ) 2 f ^ t lPel-Pe2 



In the case of bosonic neutrinos we end up with 



M b 0+ = |Mg(0 



3(£^2 — -^l) 2 + (-£ e 2 — -^el) 2 

48m 2 e 

9{E u2 - E ul f - (E e2 - E el f p el ■ p e2 



^ = + (27) 

The Gamow- Teller matrix elements are given by 

GTV 7 ./^ ^ f — E7. _l_ R n 0\r v 7 



(r = 1,2,3). 

The full decay probabilities in pure bosonic Wb and pure fermionic Wf cases can be 
written as 

W)(0 + ) = |M«(0 + )| 2 4 5D (0 + ), 

W)(2+) = \M^(2+)\% SD (2+) (29) 

and 

w b (r) = \M% ) T (r)\*i{ ISD (r) ) (so) 

where the phase space integrals are given by 

Zhsd(J w ) = — if / Jj'tt {E e i, E e2 , E u i, E u2 )F (Zf, E e i)p e iE e idE e i x 

pEi—Ef—E el pEi—Ef—E e i—E e 2 

/ F (Z f ,p e2 )p e2 E e2 dE e2 / El 2 E 2 ul dE ul (31) 

with £"^2 = Ei — Ef — E e i — E e2 — E v \ and 



2 



fj*(E e i, E e2 , E v i, E u2 ) — 



(E e \ — E e2 ) 2 [E v \ — E v2 ) 
2m% 

3(E u2 — E v i) 2 + (E e2 — E e i)' 



48m 



(Eei — E, 



c2 



\2 



Ami 



(r 


= o + ), 


(.r 


= 2 + ), 


(.r 


= o + ), 


(.F 


= 2 + ). 



(32) 
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The 2i//3/3-decay half-life is 



T ^ jr) = H^y- (33) 



3.2 Single state dominance 

The single state dominance hypothesis assumes that the 2i//3/3-decays with 1 + ground 
state of the intermediate nucleus (e.g., A=100, 116 and 128 nuclear systems) are only 
governed by the two virtual /3-transitions: i) the first one connects the ground state of 
the initial nucleus with 1+ intermediate state; ii) the second one proceeds from if state 
to the final ground state. In this case we find 



Mq+ = \M g . s X0 + )\ 2 m 2 e 



^ K f,b K f,B + L f,b L f,b + K f,b L f, b) 

3 



±(2K f > b K f > b + 2L^ b + 5K f > b L f ' b ) ?§^- 
9 E el E e2 

M% = m\ |M, S .(2+)| 2 (K*> b - l/' 6 ) 2 (l + (34) 
with K f ' b = K f m b =1 , L f ' b = L f m b =1 and 

M g . s .(r) = ±=< ^11^^111+ >< lt\\J2 T >kW°i > ■ ( 35 ) 



The value of the matrix element M g _ s .( J"") can be determined in a model independent 
way from the single /?-decay and electron capture measurements. From the experimental 
values of log ft u for the electron capture and the single (3 decay of the ground state of 
the intermediate nucleus with J 71 " = 1 + we obtain 



<iniE T * +ff *ii°* +> i = t 



3D 



9a V ftEC ' 



+ 1 1 , , 1 3D 

\<mE^Mt>\ - -fa. 

Here D = G*g 4 A /(8n 7 ). 

Within the SSD approach for the full decay probabilities we find 

w Lb (r) = \M g . s XJ n )\ 2 i f s b SD (J n )i (37) 

where 



Ei—Ef—m e 



rf.h , _ 202;/ / ^f.i 

.11 



^ssd{J w ) = — i7 / gjn(E e i, E e2} E ul , E l/2 )F (Zf, E e i)p e iE e idE e i 

"Ei—Ef—E e i r-Ei—Ef—E e i—E e 2 



X 



/ 



F (Z f ,p e2 )p e2 E e2 dE e2 / El 2 El x dE vl (38) 

Jo 



1 Because of wide range of /3-lifetimes, transitions are classified by log 10 ft values (see e.g. [E]). t and 
/ denote the measured half-life and the Fermi integral, respectively. 



9 



with 



g +{E e i, E e2 , E u i, E u2 ) 



g 2 +(E e i, E e2 , E u i, E, 



vl) 



-(K f > b K f > b + L^ b L^ b + K f ' b L f ' b ) 
3 

{Kf< b - L^f . 



(39) 



4 Characteristics of double beta decays 

In what follows we calculate the characteristics of the double beta decay mainly for two 
nuclei 100 Mo and 76 Ge for which the highest number of events has been collected in 
experiment (see Ref. [TB] and [T7j respectively). 

4.1 Double beta decay of 100 Mo 

The NEMO-3 collaboration has detected about 219 000 (0+ -> 0+)-decays of 100 Mo PS]. 
The signal to background ratio is very high S/B = 44 and the background is at the level 
of 2.5% only. All parameters of the decay: the sum of the electron energies, the energy 
of each electron and the angular distribution (angular correlation of electrons) have been 
measured. 

In the case of 100 Mo the decay proceeds mainly through the 1 + intermediate nucleus 
and the single state dominance (SSD) hypothesis should give a good approximation. 
This is also confirmed by spectra measurements in NEMO-3 experiment [TS], HH] . Since 
E m — Ei ~ Ei — Ef, the lepton energies are important in the energy-denominators (Q, 
and consequently, in the rates. 

In the SSD approximation one can calculate the probability (NME) using existing 
experimental data for the beta-decay and the electron capture of 100 Tc which is the 
intermediate dominating state. Accuracy of this "phenomenological" calculation is about 
50%, mainly because of poor experimental accuracy for the electron capture process. 

Using the SSD approximation we calculated the 2i//?/3-decay half-life of 100 Mo to 
ground state for fermionic [2] and bosonic neutrinos (see sect. 3) 

T[/M.s) = 6 - 8 10 18 years, T b /2 (0+J = 8.9 10 19 years, (40) 

so that the ratio of probabilities equals 

r (0+J = 0.076. (41) 

The ratio r Q (0^ s ) determines the weight with which the bosonic component enters the 
total rate and differential distribution [see Eq. (jl5j . For small r , a substantial modification 
of the distribution is expected for sin 2 x being close to 1. 

The higher intermediate levels can give some (basically unknown) contribution and this 
produces a systematic error in our analysis. To evaluate effect of the higher states, one can 
consider the extreme case described by the higher states dominance (HSD) approximation, 
which allows one to factorize the nuclear matrix element and integration over the phase 
space of outgoing leptons. In this case the main contribution to the 2z//5/5-decay matrix 
element comes from the transition through higher energy states (including the region 
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T [MeV] 

Figure 1: The differential decay rates normalized to the total decay rate vs. the sum 
of the kinetic energy of outgoing electrons T for 2z//3/3-decay of 100 Mo to the ground 
state of final nucleus. The results are presented for the cases of pure fermionic and 
pure bosonic neutrinos. The calculations have been performed within the single-state 
dominance hypothesis (SSD) and with the assumption of dominance of higher lying states 
(HSD). 



of the Gamow- Teller resonance) of the intermediate nucleus. Thus, the lepton energies 
in the denominators (jUJ) can be neglected (or approximated by (Ef — E i )/2 ) due to a 
large value of E n — Ei. The fermionic and bosonic 2i//5/3-decay rates are associated with 
different nuclear matrix elements [see Eq. ( 1261) and ( 1271) ]. They can be evaluated within 
an appropriate nuclear model like Quasiparticle Random Phase approximation (QRPA) 
or Nuclear Shell Model (NSM). Then, the evaluated values of 2i//3/3-decay half-life and 
ratio ^o(0^" s ) are model dependent. Contrary, the normalized differential characteristics 
are model independent for cases of pure fermionic and bosonic neutrinos. 

The energy spectra of electrons calculated in the SSD and HSD approximations are 
presented in the figs. (JTJ) and (j2J). The SSD approximation gives slightly wider spectra of 
two electrons both for fermionic and bosonic neutrinos. The spectra for bosonic neutrinos 
are softer in both approximations. In particular, the maxima of SSD and HSD spectra 
are shifted to low energies for bosonic neutrinos by about 15 % with respect to fermionic- 
neutrino spectra. This shift does not depend on the approximation and therefore can be 
considered as the solid signature of bosonic neutrino. Also the energy spectrum for single 
electron becomes softer in the bosonic case (Fig. [5]). 
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Figure 2: The single electron differential decay rate normalized to the total decay rate vs. 
the electron energy for 2i//3/3-decay of 100 Mo to the ground state of final nucleus. E and 
m e represent the energy and mass of the electron, respectively. The results are presented 
for the cases of pure fermionic and pure bosonic neutrinos. The conventions are the same 
as in Fig. CD 



In Fig. [3] we show the energy spectra of two electrons for different values of the bosonic- 
fraction sin 2 x- With increase of sin 2 \ the spectra shift to smaller energies. Due to 
smallness of tq substantial shift occurs only when sin 2 x is close to 1.0 

In Fig. H] we show the energy spectra of single electrons for different values of sin 2 x- A 
substantial change occurs at very low energies, with E kin = 0.3 MeV being a fixed point. 
For Ek in < 0.3 MeV the distribution increase with sin 2 x, whereas for Ek in = 0.3 — 1.4 
MeV it decreases. 

As we mentioned before, the rates of transitions to first excited 2^ state are affected 
by the presence of bosonic neutrino component in the opposite (to + ) way. Furthermore, 
in the SSD approximation the ratio of decay rates to the excited 2 + state and to the 0+ s 
ground state does not depend on the log ftsc value, which is not measured accurately 
enough. For the 2^/?/3-decay of 100 Mo within the SSD approximation we obtain 

^1/2 (^l") = 1-7 10 23 years (fermionic u) 

= 2.4 1 22 years (bosonic v). (42) 
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T [MeV] 

Figure 3: The differential decay rates normalized to the total decay rate vs. the sum of 
the kinetic energy of outgoing electrons T for 2i//3/3-decay of 100 Mo to the ground state of 
final nucleus. The results are presented for different values of the squared admixture sin 2 \ 
of the bosonic component. The spectra have been calculated in the SSD approximation. 



Then the ratio of the bosonic and fermionic half-lives equals 

r (2+) = 7.1. (43) 

The bosonic rate is larger in agreement with our qualitative consideration in sect. 2. 

The best lower bound on the 2z//3/3-decay half-life to excited 2f state is 1.6 10 21 
years [2D]- The current limit of NEMO-3 experiment is 1.1 10 21 years [2T] (for 1 year of 
measurements). After 5 years of measurements with the present low-radon background 
conditions sensitivity will increase up to ~ 10 22 years thus approaching the prediction in 
the case of bosonic neutrinos. Due to the large value of r even a small fraction of bosonic 
neutrinos can produce significant distortion of the standard (fermionic) spectra. 

Modifications of the spectra are opposite for the decay of 100 Mo into 2 + excited state: 
the spectra become harder with increase of sin 2 \ ( see Fig- Eland [6]). This is apparently 
related to the change of the spin of the nuclei. In the case of + — 2± transition the 
leptonic system should take spin 2 and therefore due to polarization of leptons (determined 
by V - A character of interactions) both electrons move preferably in the same direction 
(hemisphere) and two antineutrinos in the opposite direction with the corresponding Pauli 
blocking factor. In the case of bosonic neutrinos the Pauli blocking effect is reduced 
and therefore the electrons can be more aligned and consequently have higher energies. 
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Figure 4: The single electron differential decay rate normalized to the total decay rate 
vs. the electron energy for 2z//3/3-decay of 100 Mo to the ground state of final nucleus. The 
results are presented for different values of the squared admixture sin 2 x of the bosonic 
component. The spectra have been calculated in the SSD approximation. The conventions 
are the same as in Fig. [2j 



Correspondingly the spectrum becomes harder. In the case of + — + transition the total 
leptonic momentum is zero, so that the electrons move in the opposite directions. 

According to Fig. [5] even 10 % of "bosonic" admixture gives substantial distortion 
effect and this fact can be used in the future experiments. 

The angular distribution of outgoing electrons [13] can be written as 

a cos 9 2 

where 9 is the angle between two electrons. For + — + transition and fermionic neutrinos 
in the SSD approximation 

k (0+ s ) = —0.627 (fermionic neutrino). (45) 

(The HSD approximation gives similar number: —0.646.) Notice that the preferable 
direction is 9 = 180° when electrons move in the opposite directions. The configuration 
with the same direction of two electrons is suppressed. For bosonic neutrinos we find 

k\0+ s ) = -0.344 (bosonic neutrino). (46) 
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1.2 




T [MeV] 

Figure 5: The differential decay rates normalized to the total decay rate vs. the sum 
of the kinetic energy of outgoing electrons T for 2z//3/3-decay of 100 Mo to the excited 
2^ state of final nucleus. The results are presented for different values of the squared 
admixture sin 2 \ of the bosonic component. The spectra have been calculated in the SSD 
approximation. 



(The HSD approximation gives —0.422.) So, the configuration with the same direction 
of electrons is less suppressed and the distribution is more isotropic (flatter) than in the 
fermionic case. 

4.2 76 Ge double beta decay 

The statistics of 76 Ge decays is about 113000 events, the background is rather high, S/B 
= 1.3, and only the sum of two electron energies is measured [TTj . The systematic error can 
be as large as 10% and the main source of the error is the background. One has to estimate 
this background independently and make subtraction. So, one can shift the spectrum and 
its maximum within the error. Furthermore, the energy spectrum of two electrons starts 
to dominate over the background above 0.7 MeV which means that the maximum of the 
spectrum is not observed. The advantage of 76 Ge is that there is practically no difference 
between the results of HSD and SSD approximations for the energy distributions because 
the nearest 1^ state of the intermediate nucleus is lying high enough. Thus, one does 
not need to make assumptions about SSD or HSD. In this way the conclusion does not 
depend on the nuclear structure details. 
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Figure 6: The single electron differential decay rate normalized to the total decay rate vs. 
the electron energy for 2i//?/3-decay of 100 Mo to the excited 2f state of final nucleus. The 
results are presented for different values of the squared admixture sin 2 x of the bosonic 
component. The spectra have been calculated in the SSD approximation. The conventions 
are the same as in Fig. [2j 



In the HSD approximation, evaluating the phase space integrals and nuclear matrix 
elements within the proton-neutron QRPA we find 

r (0+J = 0.0014. (47) 

This smallness is related to a large extend to high energies of the intermediate states, 
E m — Ei in comparison with leptonic energies restricted by the energy release E\ < (Ei — 
Ef)/2: Ei <C E m — Ei. According to ([9]) the factors K h m , L b m and consequently the rate 
are zero in the limit Ei = 0. In the lowest approximation we obtain 

K b T b ^ Kjvg ~ E vl ) ± [E e2 - geijj (dH , 

{E m ~ Ei) 2 

(where plus sign is for ^-factors). Then the ratio of the rates can be estimated as 



where e\ is the average energy of the lepton. Taking parameters of the 76 Ge -decay we 
find ro ~ 10 -3 in a good agreement with the calculations in QRPA. 
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T [MeV] 

Figure 7: The differential decay rates normalized to the total decay rate vs. the sum of 
the kinetic energy of outgoing electrons T for 2z//3/3-decay of 76 Ge to the ground state of 
final nucleus. The results are presented for the cases of pure fermionic and pure bosonic 
neutrinos. The calculations have been performed with the HSD assumption. 



In Fig. [7] we show the normalized distributions of the total energy of two electrons for 
pure fermionic and bosonic neutrinos. As in the case of 100 Mo, the decay with bosonic 
neutrinos has softer spectrum. The energy distribution of single electron is shown in Fig. [H] 

Due to a small value of r (0+ s ) a substantial effect of the bosonic component should 
show up only for sin 2 \ being very close to 1: (1 — sin 2 x) 2 ~ 10r (0+ s _). So studies of the 
spectra are not sensitive to sin 2 x- In contrast, the total rate of the 76 Ge decay gives a 
strong bound on sin 2 \- 

5 Bounds on bosonic neutrinos 

One can search for/restrict the bosonic or partly bosonic neutrino using total rates, ratios 
of rates of the transitions to the excited and ground states, energy spectra, and angular 
distributions. Let us evaluate the bounds on sin 2 x that can be obtained from the existing 
data using these methods. 

As follows from our general discussion in sec. 3, for + — > + transitions: r <C 1. For 
nuclei with small r the best bound on bosonic neutrino fraction can be obtained from 
the total rates. A modification of the spectrum due to presence of bosonic component is 
small. In contrast, the strongest modification of the spectrum is expected for the nuclei 
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Figure 8: The single electron differential decay rate normalized to the total decay rate 
vs. the electron energy for 2i//?/?-decay of 76 Ge to the ground state of final nucleus. E and 
m e represent the energy and mass of the electron, respectively. The results are presented 
for the cases of pure fermionic and pure bosonic neutrinos. The calculations have been 
performed with the HSD assumption. The conventions are the same as in Fig. [2j 



with large tq. This is true, e.g., for + — > 2 + transition, where ro ^> 1. 



1) Method 1: Comparison of the predicted and measured half- life times. Using (I15p 
we can write 



sin 2 x 



1 



l + r 



1 - 



T V 2 ! I 

=5 - r 1 



\ rpexp 
\ J V2 



1 1/2 
rpexp 
1 l/2 



(50) 



where r = T^ 2 /T^ 2 , T(, 2 (T^ 2 ) are the theoretically predicted life-times for fermionic 
(bosonic) neutrinos and Ty 2 is the experimentally measured life-time. In the case of 
agreement between the measured and the predicted (for fermionic neutrinos) life-times, 
we can use ( 15*01 to establish the bound on parameter sin 2 x- 



sin x < 



l + r 



1 - 



\ 



f—min 



1/2 



rpexp— max 
1 l/2 



-r 1 



rp f—min 

1 1/2 
rpexp— max 
1 l/2 



(51) 



Here, T^ 2 mm and Ty 2 max are, respectively, minimal theoretical value within a considered 
nuclear model (e.g., QRPA and its modification, NSM) and maximal experimental value 
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of the permitted experimental range of the 2u (3 /3-decay half-life. For r C 1 and r smaller 
the relative accuracy of determination of T^ 2 /T^J 2 the terms proportional to r in (l5*Tj) 

can be omitted. Then we get sin 2 x < (1 — J T[j 2 mm j r p^jp- max y 

Apparently, this method requires knowledge of the nuclear matrix element, and as 
we mentioned above, reliable estimations can be done for some nuclei e.g., 100 Mo and 
116 Cd assuming SSD hypothesis. For some other nuclear systems nuclear models have 
to be considered. The two basic approaches used so far for the evaluation of the dou- 
ble beta decay matrix elements are the QRPA and the NSM. For the 2^/3/3-decay of 
76 Ge the predicted half-lives are 7.7 10 20 - 1.4 10 21 years (QRPA) [22] and 1.15 10 21 
years (NSM) [23]. The experimental half-life (average half-life value is (1.5 ± 0.1) 10 21 
years [21]) is in rather good agreement with the theoretical ones for fermionic neutrino 
within uncertainty characterized by the factor ~ 2 (see [22])- For pure bosonic neutrinos 
r o(0^ s .) ~ 10~ 3 (QRPA) and therefore for the half-life time we would have T\, 2 « 1.5 10 24 
years, which is in contradiction with the experimental value. So, purely bosonic neutrino 
is certainly excluded. 

The axial-vector coupling constant qa is a significant source of uncertainty in the 
theoretical calculation of the 2z//9/?-decay rate, which is proportional to g\. The commonly 
adopted values are Qa = 1-0 (by assuming quenching in nuclear medium) and qa = 1-25 
(as for free nucleon). This gives about 1.5 uncertainty in NME's. 

For factor 2 uncertainty in NME we obtain factor 4 uncertainty in Ty 2 . Therefore 

taking T[, 2 ~ Ty 2 , we can put the bound 

rpf—min 

llll > I (521 

rpexp—max a' V / 

J l/2 

Then, eq. ( 15*11 gives 

sin 2 x< 0.50. (53) 

Notice that uncertainty in T[, 2 (and not r ) dominates in this bound. 

We can also use the half-life time of 100 Mo. Here ro(0+ s ) is much larger ( |4"T1) but the 
accuracy of calculations of NME is better. Taking SSD approximation we can calculate 
the half life with 50% accuracy: Ty 2 = (6.84 ± 3.42) 10 18 years [13]. This value is in 
agreement with NEMO-3 value, T*T = (7.11 ± 0.54) 10 18 years [16J. Plugging these 
numbers into (I5*T1) we obtain for r (0+ s ) = 0.086 

sin 2 x < 0.34. (54) 

Notice that the accuracy of predicted half-life value is connected with experimental ac- 
curacy for EC (electron capture) half-life of 100 Tc [25]. This accuracy can be improved 
in the future experimental down to ~ 10% and correspondingly, the sensitivity to sin 2 x 



2 In ref. [25] Mo enriched to 97.4% was used and the main background was connected with X-rays from 
different Tc isotopes which were produced in the sample due to (p,n) and (p,a) reactions on different 
Mo isotopes, from 92 Mo to 98 Mo; see Table II in [25 . If one uses Mo enriched to 99% (or more) then 
the mentioned above background would be much lower and the accuracy of the measurement would be 
several times better. 
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can reach ~ 0.1. Unfortunately, there is only one (not very precise) EC measurement for 
100 Tc and thus the above limit on sin 2 x is not reliable enough. 

Even stronger bound can be obtained from studies of 116 Cd -decay. Recently a precise 
estimation of half-life value based on the SSD approximation and information from the 
116 Cd(p,n) reaction was obtained: T/ /2 = (2.76 ±0.12) 10 19 years [26]. This prediction 
is in a very good agreement with experimental value (The experimental average is (2.8 ± 
0.2) 10 19 years [24J). Using these results we obtain from (JBTl) 

sin 2 x< 0.06. (55) 

It should be noticed that the result of ref. [26] substantially differs from the earlier es- 
timation T/ /2 = (1.1 ± 0.3) 10 19 years [H] (also based on SSD and measured value of 
electron capture rate of 116 In [2Z])- This result disagrees with the experimental value and 
could be interpreted as the effect of partly bosonic neutrino with sin 2 x ~ 0.4. 



2) Method 2: Measurements of the differential characteristics of the decays: shapes 
of the energy spectra (sum energy and single electron energy) and angular distribution. 
Such information is provided now by NEMO-3 for 100 Mo, 82 Se, 116 Cd, 150 Nd, 96 Zr and 
48 Ca. In the future the results for 130 Te will be also available pU M, ES EE]. In this 
method one compares the experimental and theoretical energy spectra as well as the 
angular distribution. In practice one should perform the statistical fit of the spectra by a 
general distribution (|T6|) with sin 2 x being a free parameter. As we have seen the spectral 
method has substantial sensitivity to sin 2 x for nuclei and transitions with large r . That 
includes 100 Mo, as well as transitions to the excited states. 76 Ge with very small r has 
no high sensitivity. 

a) Let us consider first the energy spectra of 0+ s — > 0+ s decay of 100 Mo [16]. In the 
present paper we will not perform detailed statistical analysis of the spectra, postponing 
this to the time when measurements will be finished and all careful calibrations will be 
done. Instead, we give some qualitative estimates. There is a reasonable agreement 
between the predicted energy spectrum of two electrons and the experimental points. 
Therefore we can certainly exclude the pure bosonic case (sin 2 % = 1). Furthermore, 
comparing the results of Fig. [3] (essentially, the relative shift of the maximum of spectrum) 
with the experimental spectrum we can put the conservative bound sin 2 ^ < 0.6. In fact, 
there is no ideal agreement between data and theoretical spectrum. A better fit can be 
obtained for sin 2 x ~ 0.4 — 0.5. 

b) Let us comment on the single-electron energy spectrum from 100 Mo decay. The data 
reasonably well agree with the predictions from the fermionic SSD mechanism, but some 
difference exists between the data and the fermionic HSD-mechanism predictions. From 
this it was concluded that the SSD mechanism is more relevant here [18, 19J. Comparing 
the experimental data and spectra for partly bosonic neutrinos (Fig. H]) we obtain: sin 2 x < 
0.7. 

Notice that the SSD spectrum does not show an ideal agreement with data either. 
There is some discrepancy, especially in the low energy region (E = 0.2 — 0.4 MeV). That 
could be explained by the effect of partly bosonic neutrinos with sin 2 x ~ 0.5 - 0.6. 
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Complete analysis of all existing NEMO-3 information (energy and angular distribu- 
tions) using e.g. maximal likelihood methods, will have a higher sensitivity to sin 2 %. 
However, it is difficult to expect a better bound than sin 2 x ~ 0.4 — 0.5, mainly because 
of the existing disagreement between the data and Monte Carlo (MC) simulations. In 
fact, it can be just some systematic effect connected to the present poor understanding of 
response function of the detector. If in future the NEMO experimental data turn out to 
be in much better agreement with the MC-simulated spectrum, the sensitivity to partly 
bosonic neutrino will be improved down to sin 2 x — 0-2 — 0.3. 

3) Method 3: Determination of the ratios of half-lives to excited and ground state, 



1 l/2\ J 

-^1/2(^5.) 



separately for fermionic and bosonic neutrinos. For 2t//9/?-decay of 100 Mo the ratio can be 
calculated rather reliably using the SSD- approximation. The advantage of this quantity 
is that the EC amplitude, [(A,Z) — > (A,Z+1) transition], which is not well determined, 
cancels in the ratio ( |56l) . 

For 100 Mo the transitions to the ground 0+ s and excited Of states were detected, and 
in fact, some discrepancy has been observed. The corresponding experimental ratio r* 
equals 

r e V(0+)~ 80 (57) 
(NEMO-3 results (El EI]), whereas within the SSD approach the calculated ones are 

r*(0f) ~ 61 (fermionic v) 

~ 73 (bosonic v). (58) 

A bosonic neutrino fits the data slightly better but the differences are probably beyond 
the accuracy of the SSD assumption. Still it is also necessary to improve statistics in 
measurements of the transition to excited 0i state. 

Contrary to the case of + excited state, the ratio of 2^/9/?-decay half-lives to excited 2 + 
and ground state is expected to be strongly different for bosonic and fermionic neutrinos. 
Using the SSD approximation for the 2i//9/9-decay of 100 Mo we found 

r*(2f ) ~ 2.5 10 4 (fermionic v) 

~ 2.7 10 2 (bosonic v). (59) 

The 2i//3/3-decay of 100 Mo to excited 2f state has been not measured yet. Using the best 
experimental limit on the half-life found in [20] we get 

r e * xp (2+) > 2.2 10 2 . (60) 

This bound is close to the bosonic prediction. A further experimental progress in mea- 
suring this nuclear transition will allow one to analyze also the case of partially bosonic 
neutrino, and therefore is highly required. 
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6 Conclusions 



A study of the double beta decay can provide a sensitive test of the Pauli exclusion prin- 
ciple and statistics of neutrinos. (Notice, that relation between the statistics of neutrinos 
and possible (small) violation of the Pauli principle is an open issue.) Appearance of 
the bosonic component in the neutrino states changes substantially the total rates of the 
decays as well as the energy and angular distributions. We find, in particular, that the 
ratio tq(0^ s ) of the rates to ground state for bosonic and fermionic neutrinos, is < 1CT 3 
for 76 Ge and 0.076 for 100 Mo, which excludes pure bosonic neutrinos. For transitions to 
2 + excited states tq(2 + ) ^> 1, in particular ?"o(2^) ~ 7. However, this 2i//3/3-decay channel 
has been not measured yet. 

We have introduced phenomenological parameter sin 2 \ that describes the mixed statis- 
tics case of partly bosonic neutrinos. The dependence of the energy spectra and angular 
correlation of electrons on sin 2 x has been studied. The bound on sin 2 x can be obtained 
by comparison of the predicted and measured total rates of the decays. In spite of the 
big difference of the rates for fermionic and bosonic neutrinos, this method does not give 
strong and very reliable bound on sin 2 x due to uncertainties in NME's. The conserva- 
tive upper bound sin 2 x < 0.5 is found using the 100 Mo and 76 Ge results. Much stronger 
bound, sin 2 x < 0.06, is obtained from recent studies of 116 Cd, however this bound re- 
quires further checks. 

The method based on the study of the normalized energy and angular spectra is less af- 
fected by uncertainties in the NME's. The transitions with large ^(J 71 ^) have the highest 
sensitivity to spectrum distortions and therefore sin 2 x- Using the data on the 0+ s — > 0+ s 
transition of 100 Mo we obtain the bound sin 2 x < 0.6. In the future this bound can be 
improved down to sin 2 x ~ 0.2. The 0+ s — > 2f transition with 7"o(2+) ~ 7 can give much 
stronger bound, but here new, more sensitive experimental results are needed. We find 
that modification of the energy spectra due the presence of the bosonic components is op- 
posite for 0+ s — > 0+ s and 0+ s — > l\ transitions: for 0+ s — > 0+ s the bosonic component 
leads to softer spectrum whereas for 0+ s — > l\ transitions to harder spectrum of elec- 
trons. Also the presence of bosonic component leads to flatter angular (cos 9) distribution. 

Strong bound (potentially down to sin 2 % ~ 0.1 — 0.05) might be obtained from mea- 
surements of ratios of the decay rates to the 2+ excited and ground state. However, this 
requires further experimental progress. 

We note that currently there are no restrictions on the admixture of bosonic component 
from the BBN. However, as it was indicated in [7] the future BBN studies will be able to 
constrain the fermi-bose parameter to k > 0.5. The bound on parameter sin 2 x < 0.6 
from the 2z//3/3-decay results in k > —0.2. 

In conclusion, the present data allow to put the conservative upper bound on the admix- 
ture of the bosonic component sin 2 x < 0.6. With the presently operating experiments 
this bound might be improved down to 0.2. In future one order of magnitude improvement 
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seems feasible. 
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